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Abstract 

We extend our study of thermodynamics of a Kubo particle to temperatures smaller 
than the mean interlevel spacing. We obtain the distribution functions of spin susceptibility 
and heat capacity for Poisson and Wigner-Dyson level statistics. We evaluate the line shape 
of the Knight shift due to spin effects both in a single particle and for the ensemble average 
and compare it with orbital and spin-orbit contributions. 
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1 . The spin susceptibility of an ensemble of small particles with Poisson level statistics 
was first evaluated by Kubo [1]. Building on the absence of level correlations, he found the 
exact result both for the odd and even electron particles. Gor'kov and Eliashberg argued [2] 
that disordered particles and particles with diffusive boundaries should have Wigner-Dyson 
(WD) level statistics [3]. Denton, Miihschlegel, and Scalapino [4] have noticed that Kubo's 
results follow from a two- and three-level model in the even and odd cases respectively and 
evaluated the electron spin susceptibility and heat capacity for the WD circumstance. These 
developments have been extensively reviewed [5]. 

2 . Despite the realization of the statistical nature of spin response, the authors of these 
works have been concerned only with the average susceptibility. On the other hand, the test 
of their predictions relied mainly on the measurements of the Knight shift [6]. Obviously, 
the latter should be broadened due to the distribution of spin susceptibilities. Although the 
inhomogeneous broadening has been of concern [5], it was attributed largely to the particle 
size distribution. Below we obtain the distribution of the spin susceptibility using the model 
of Ref. [4] and discuss its effect on the inhomogeneous broadening of the Knight shift vis- 
a-vis the fluctuations of the wave function in the particle [7] . We compare this effect to the 
inhomogeneous broadening due to orbital and spin-orbit effects [8]. We also evaluate the 
distribution function of the heat capacity. 

3 . In what follows, we set the Curie susceptibility, Xc = f^lP ■> where (5 = {k B T} 1 and 
jilg is the Bohr magneton, to unity. In these units, the Pauli susceptibility is % P = £ , where 
r 1 = /3 8 » 1 and 8 is the mean interlevel spacing. For even electron particles, the partition 
function of the two-level system - the Fermi level, taken at zero, and the nearest unoccupied 
level at a distance A - in the magnetic field is given by [4], [5], 

Z = 1 + 2( 1 + cosh \2Pfi B HUexp [- fiAj + exp (- 2/Jzl) , (1) 
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The distribution function for the susceptibility is obtained via integration of the Dirac delta- 
function whose argument is obtained from Eq. (1), 



p(x) = 



dAd 



[z- TT ^ T? )/'M, * = exp(-/fc), (2) 



where P (zl) is the nearest level distribution function given by [3], [5], 
P (A) = iexp I- j] , Poisson, (3a) 



P(4) = f Jrexp|-|^gjj, orthogonal, (3b) 
for the two cases considered here. 

4 . The integral (2) can be trivially evaluated with the help of the formula, 

d(f{x]) = Z\f(X,)\ l d(x-X}), (4) 

where X t are the roots of the equation / (X,) = . The result of exact integration, however, 

is quite messy as well as partially superfluous. A simplified expression can be obtained on 
the basis of the following argument. The mean susceptibilities are (%) ~ e and (y\ ~ e 2 in 

the Poisson and orthogonal cases respectively. Using the result of exact integration, it can 

r- r max 

be also shown that d%P \X) is ~ £ and ~ e 2 respectively, while d%p {%) ~ 1 in 

J{x) ' Jo 

either case, where # max ~ 1 is the maximum value of the susceptibility found from Eqs. (2) 
and (4). We further notice that the accuracy of the distributions obtained from the two-level 
model is reduced near % m . dX . Indeed, the n 'th moments scale as £ and e 2 regardless of n 
(see Appendix I), meaning that l^x - ^ ~ Xc~ 1 Xp m the Poisson and ~ Xc~ 2 Xp m tne 

orthogonal cases respectively which, functionally, can be understood as %c multiplied by 
the probability of finding the next level within k B T of the Fermi level. The first of these 
relationships is consistent (see Appendix II) with the formal extension of our derivation [9] 



f ( A\2 
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for ft d « 1 to (3 d» 1 . It indicates that the tails of the distributions are defined by the higher 
moments. However, as shown in Appendix I, such moments are not accurately defined by 
the two- level model, even for ft 8» 1 . 

5 . We conclude that whereas the distributions have long tails upward % ~ 1 , the two- 
level model is insufficient for their accurate description. Since also the susceptibilities larger 
than the mean have a small probability, a good approximation is achieved by taking t « 1 in 
Eq. (2). Subsequent verification confirms that the expressions thus obtained, 

P(%) = ^ £ [^x) > Poisson, (5a) 

2 I I I \\ 2 \ I \ 

p(%) = -^exp -itfjeln^jj ln||^j , orthogonal. (5b) 

coincide, within the accuracy of the model, with the result of exact integration (see Fig. 1). 
Both expressions have integrable singularities at % — > . Clearly, the small susceptibilities 
reflect the energy gap at the Fermi level. 

6 . For odd number of electrons, one must consider the three-level model and = 1 , 

irrespective of the level statistics [1], [4], [5]. From the latter fact alone one concludes that 

the distribution function is narrower than for the even electron circumstance. Formally, the 

susceptibility distribution function can be derived by averaging the delta function with the 
distribution function P [A, ZT) for the two levels [4]. We notice, however, that in the same 

approximation as above one finds the functions that diverge at ~ 1 similarly to Eqs. (5) at 

zero, with the extra factors of e and e 1 respectively (see Appendix III) that can be traced to 

the probability of finding a level within k B T below, as well as above, the Fermi level. The 

distribution tails extend to the values of the susceptibility of order of the mean. 

7 . The distribution function of the heat capacity has the same functional form for even 
and odd number of electrons and differs only in numerical coefficients. Consequently, we 



5 



limit our discussion to the even case here. From the partition function in Eq. (2) for H = 
we find for the distribution function, 



p(c) = 



dAS 



< At ^ 



c ~ W 2 T747T7 P M ' t = exp (" PA ) ' (6) 



where C is measured in units of k B . Making the same approximation as before, we have, 
dAS (y - 4exp (- pA))p (A) (6a) 



P (y) « £ 2 

Jo 



where y = Ce 2 . As a result, replacing 8 — > 4 , one finds the same expressions for p (y)e^ 2 
as for p [x] in Eqs. (5). Therefore, the heat capacity distributions have tails toward C ~ 1 , 
while the mean values are ~ e and ~ £ 2 for the Poisson and orthogonal cases respectively. 

8 . To test the predictions for the distributions of spin susceptibility and heat capacity 
experimentally, one must be able to measure the respective quantities in a single system. 
The rings and grains studied in recent experiments on orbital magnetic response [10] are a 
candidate but they are perhaps too large to achieve the condition « 1 and perform a 

reliable measurement. The ensembles of particles studied earlier [5], as well as ensembles 
of rings [10], would yield only average quantities. The Knight shift broadening, however, 
would reflect the susceptibility distribution. For this reason, we now turn to its discussion 
in the light of the above results. As is well known [11], an electron creates a magnetic field 
at the nucleus site f , 

H {f] = ^d\ ¥ {r)\ 2 , (7) 

where the over-bar denotes the quantum-mechanical average of the spin state. Performing 
now the thermodynamic averaging, we notice that in the limit (fid) ' « 1 the electrons that 

occupy the levels below the Fermi level do not contribute (for the even electron case) since 
we neglect their excitations and their spins are pair- wise opposite. The two electrons at the 
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Fermi level are allowed to jump to the next level in the two-level model. Of the six possible 
states, only two are magnetic yielding, 



9{{r) = 



S7UjLl B 


f 

exp 

I 


( -P(A-2h b h) 


-exp 


( -P(A + 2» B H^ 


3 1 + 4exp (- pA) +exp (- 20 A) 



\2 2 

¥o(r)\ + Yi(r) , (7a) 



where the wave functions are those of the Fermi level and the next level. Expanding in H , 
one finds the expression which, aside from the WFs and numerical coefficients, coincides 
with the expression for the total magnetic moment of the particle obtained from the partition 
function in Eq. (1). The same is true for odd electron particles. Therefore, we conclude that 
in a given system the Knight shift broadening should reflect the WF behavior, while on the 
average it should be a combined level statistics and WF effect. In particular, in the WD case 

I |2 

the line shape of a single particle should follow the Poisson distribution of | \ff\ [7] (notice 

that since WFs of different energy states are distributed independently [7], the contributions 
from the electrons with opposite spins « | yr (r) - \jf x (r) j average out to zero in every 

particle). On the average, however, it can be shown that the WF distribution does not effect 
significantly the result for the line shape anticipated on the basis of level statistics alone. In 
the orthogonal case, for instance, the functional dependence of the line shape should follow 
closely that of Eq. (5b) in the limit e « 1 (see Fig. 2). In classically regular systems, on the 
other hand, the WF behavior can be quite complex [7] and we will leave the discussion of 
the effect this might have on the Knight shift to future work. 

9 . We will now discuss the relationship of the above results to other mechanisms and 
regimes for the Knight shift broadening. Firstly, Efetov and Prigodin [12] have considered 
the same problem. We believe, however, that their model, wherein the chemical potential is 
presumed fixed and the level broadening is smaller than the level spacing, is unrealistic. In 
the preceding paper [9], we discussed the susceptibility fluctuations in the opposite limit, 
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ft 5 « 1 . It should be emphasized that contrary to the wide-spread belief found in the small 
particle literature [5], the discreteness of the spectrum is manifested even in this limit [9], 

* W 'l 4«0f. WD' <8) 

The analysis in Ref. [9] can be extended to higher moments (/' = d e f ), 

[z-(x) 



where (jq = Xp [9]- It can be argued that in the considered limit the largest contribution to 
Eq. (9) comes from pair-wise decoupling of the n -point level density correlation function. 
Elementary combinatorics then yields a Gaussian distribution of the susceptibility. Thus, a 
narrow Gaussian distribution centered around % P transforms to a broad distribution whose 
tail extends to ~ % c as the temperature decreases below the mean level spacing. It might be 
possible to test the susceptibility distribution for f3 S« 1 by measuring the line shape of the 
Knight shift as well. 

10 Previously, we have predicted [8] inhomogeneous broadening of the Knight shift 
whose origin is the large orbital currents in disordered metals. In the presence of spin-orbit 
scatterers, such currents can also couple to spin degrees of freedom [8]. In small particles, 
the local susceptibility fluctuations were predicted to be of the order of magnitude of the per 
unit volume Pauli (Landau) susceptibility. These effects should be compared to the Zeeman 
effects studied here. Clearly, for j58 « 1 the orbital and spin-orbit effects will be larger (the 
former should be identifiable, however, in quasi- 2D systems in an in-plane magnetic field). 
At this point, we cannot definitively predict the magnitude of orbital and spin-orbit effects 
for f3 8 » 1 since the derivation of Ref. [8] assumes a fixed chemical potential. We argued 
in Ref. [9] that such an assumption would definitely not yield considerable difference with 
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the canonical ensemble only in the opposite limit, f5 8« 1 . But since Ref. [8] predicts little 
temperature sensitivity of the fluctuations (once the length L T = \JtiD I T , where D is the 
diffusion coefficient, becomes larger than the particle size), we can extrapolate its result to 
p~ 1 ~ 8 . Notice that the upper bound on the applicability of linear response is different for 

the orbital and spin mechanisms. For the orbital response, we refer to its discussion in Ref. 
[8] where it was shown that the magnetic moment of a particle can be as large as /n B (k F C^ , 

where d is the electron mean-free-path and k F is length of the Fermi wave vector. For the 
spin response, the upper bound of the linear regime is [l B H < /T 1 . Therefore, for f38» 1 , 
the maximum magnetic moment is ~ % C H < \x B ; this estimate holds up to /T 1 ~ 8 . Clearly, 
a more meaningful comparison requires a better understanding of orbital effects in the limit 
08» 1 . 

11. In conclusion, using the two- and three-level model, we determined the distribution 
function of the spin susceptibility for even- and odd-electron particles in the limit 8 » 1 . 
In the even case, the approximate form of the distributions for the Poisson and orthogonal 
level statistics are given by Eqs. (5). They can be characterized as being roughly power-law 
functions with integrable singularities for small susceptibilities and tails extending to ~ % c 
for large susceptibilities. The main difference between the two is in the prefactors e and e 2 
respectively, where £= ' « 1 . In the odd case, the distribution functions are narrowed 
further by additional factors of e and e 2 respectively and are peaked at ~ % c . Although in 
either case the distribution function narrows as the temperature is decreased, the tails are a 
distinct feature that should be experimentally identifiable. The larger broadening in the even 
case is a reversed effect relative to the effect expected from the particle size distribution [5]. 

12. We enjoyed several illuminating discussions with B. Goodman. We used Theorist 
and Mathematica for Power Macintosh for some of the analytical calculations and plotting. 
This research was supported by the Army Research Office Grant Number DAAL03-93-G- 
0077. 
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Appendix I 



Here we evaluate the moments of the susceptibility distribution function obtained in 
the two- and three-level approximations for a system with Poisson level statistics. Based on 
their comparison, we shall argue that the few-level model is less accurate for the tails of the 
distribution. The n 'th moment (since (%) ~ e , it can be neglected) is given by, 



(r)= 



dA 



St 



1+4? + 



dt 



St 



l+4t + t 



t s- 1 



(AI-1) 



where we have used the notation of Eq. (2) and Eq. (3a). In the large n limit, the multiplier 

f~ 1 can obviously be neglected. Although we could not find a closed form for the integral, 
it is clear that it scales roughly as a" , where a > 1 and is close to ^ , the maximum value 

of the susceptibility in the two-level model which is achieved at t = 1 (A = ). This result 
is consistent with the evaluation of (x") with the distribution function (5a), which validates 

the approximation made in the derivation of the latter. Notice that in the orthogonal case one 
can use the first equality in (AI-1) and the distribution function (3b) to obtain (x") ~ £ 2 ■ 

Using the three-level partition function from Halperin's review [5], we find that the 
lowest order correction to the two-level moment is given by, 



(r\ = 







( St ] 




dAd£ 








yl +4t + t 2 J 



I 

(1 + t'JP (A, A') = (r) 2 f dA'{\ + t'JP (A') , (AI-2) 

J 



where we used P [A, zl') = P (A}P (zl') for the Poisson ensemble [1]. By means of scaling 
the integration variable and series expansion we obtain, 

(AI-3) 



{X n \ = {x%e[ dx{\ + e -*)Y~ = (r) 2 1 +e± 
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The sum in (AI-3) can be evaluated in a closed form in terms of hypergeometric functions. 

For our purposes, however, it is sufficient to notice that in the large n limit it is limited by 

Tin and 2" from below and above respectively [13]. This demonstrates that, for any e , 

there will be a large correction to a high cumulant from the account of more levels near the 

Fermi surface. Therefore, we conclude that the tails of the distribution functions are not as 

well defined in the few-level model. This conclusion can be confirmed by a simple physical 

argument. Indeed, we notice that the maximum value of the susceptibility changes with the 
account of more levels. For instance, it is ^ in the two- and j in the three-level models. 

These values are found from A = and A = A" = respectively. In other words, the tails 
are due to the systems with several levels within k B T of the Fermi level, which, in turn, 
indicates that a large number of levels must be taken into account for their characterization. 
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Appendix II 



Here we detail a formal extension of our previous derivation [9], valid for j5 8 « 1 , 
to f$8» 1 . We also show that it is equivalent to the derivation of Ref. [12]. We begin with 
the expression for the n 'th moment of the susceptibility, 



X 



^-Zcosh-M 



,2n 



(AIM) 



where we have used the diagonal approximation which gives the largest contribution in the 
limit f3 8 » 1 . Converting (Al) to an integral on energy with the level density 1 , we find, 
upon some manipulations and integration [13], 



In Q" - 1 



i r \ = 4^— r dxcosh- 2n (x) = -4 

N 1 2 2n - l d J-~ { 1 2"- 1 



}■ 



(2n-l)\ 



_ - 1 

, AC AP 



(AII-2) 



To obtain the local moments, one must also Poisson-average the n 'th power of the squared 
WF, which yields an extra factor of n ! . The latter is illustrated for the model considered in 
Ref. [12] wherein /T 1 = but a finite level broadening is assumed. From Eq. (4.7) in the 
second of Refs. [12], we obtain in the diagonal approximation and with the help of simple 
integration [13], 



r = 



..In 



de 



dv 



Jo 



yl 2k\v 



e 2 + (7/ 2) 



exp - v 



I 2 \ n 



(„-!)! 



I 1 ^ 



v 2a ; 



n-m \2n - 2 ! 



(AII-3) 



(n-l)l 



which coincides with Eq. (4.21). Moreover, the characteristic function formalism allows to 
restore the distribution function from the moments [14]. For instance, a trivial calculation 
using (AII-3) yields the following distribution function: 
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(AII-4) 



where 9 (x) is the step-function and x = % I % P . The second term in Eq. (AII-4) obviously 

reproduces the limiting results given by Eqs. (4.16) and (4.18) in the second of Refs. [12], 
while its first term coincides with the limit a — > of Eq. (4.6). In other words, once the 
WF distribution is known [7], the derivation becomes transparent and the application of the 
non-linear a -model is unnecessary, at least for x> a . More importantly, we stress the 
formal nature of the above derivation in the limit f3 8» 1 since the first of equalities in Eq. 
(AII-1) is valid only in the opposite limit f3 8 « 1 . Furthermore, its results, given by Eqs. 
(AII-2) and (AII-3), are clearly statistics-independent, while the correct results (see text and 
Appendix I) clearly are. It helps, on the other hand, to establish a connection between our 
results for 8 « 1 [9] and the results of a few-level model which is appropriate for f3 8» 1 
and is a subject of this work. 
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Appendix III 

Here we derive the susceptibility distribution function for the odd-electron particles 
with Poisson level statistics. Using the partition function in Halperin's review [5], we find, 




(AIII-1) 



where we have used Eq. (3a) and the approximation t, f « 1 . Without the last term inside 
the delta-function, one recovers p = 8yx - 1 ) , the Curie susceptibility. The structure of 
the second integral is as follows. Without the account for the [t + fj term, one would have 
£ 2 i^% - 1 ) / 8 j for x > 1 > but the integral on second variable would yield a logarithmic 

divergence at zero. The account of this term, however, shows that the divergence is cut-off 
by [%- 1J/8 . Consequently, we find, 

P (X) = ~ £ 2 ((X ~ 1 ) / 8)" 1 In (( X - 1 ) / 8) (AIII-2) 

which, up to the extra factor of e , has virtually the same behavior at 1 as the distribution 

function (5a) at zero. Notice that in the three-level approximation used here, the maximum 
susceptibility is §■ . 
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Figure Captions 

Figure 1: p(%) given by Eq. (5b) (solid line) and p y%) obtained as a result of exact 
integration in Eq. (2) (dashed line) for e = 0.2 . 

Figure 2: p [x) given by Eq. (5b) (solid line) and Jvexp (- v)v~ l p [y] (dashed 
line) for £=0.2 . 
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